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ABSTRACT 

We report the discovery of a narrow-line QSO (NLQSO) RX J1334.2+3759 with a 
steep soft X-ray spectrum. Soft X-ray emission from the NLQSO is highly variable. 
Changes in the intensity by a factor of ~ 2 have been detected in the ROSAT PSPC 
observations ofRX J1334. 2+3759 on time scales of ~ 20000-40000 s. Rapid variability 
' events have also been observed from RX J1334. 2+3759. The most extreme variable 

c"^i \ event has — (1.95 + 1.02) x 10 42 erg s -2 corresponding to a change in intensity by a 

£NJ . factor of ~ 4 within ^ 400 s. The PSPC spectra of the NLQSO can be well represented 

by a power-law of photon index, T x ~ 3.8 modified by an absorbing column local to the 
source (ANh ~ 3.3 x 10 20 cm -2 ) over and above that due to our own Galaxy (Nh = 
7.9 x 10 19 cm -2 ). The intrinsic soft X-ray luminosity of RX J1334. 2+3759 is estimated 
to be ~ 2.2 x 10 44 erg s _1 in the energy band of 0.1 — 2.0 keV. The optical spectrum of 
r"| . RX J1334. 2+3759 is typical of the NLS1 galaxies showing strong Balmcr H/3, Ha, and 

O ^ forbidden line of [O IIIJA5007. Fe II multiplets, usually present in the optical spectra 

of NLS1 galaxies, are also detected from RX J1334. 2+3759. Decomposition of the 
H/3 and Ha line profiles shows the presence of narrow (FWHM ~ 880 km s -1 ) and 
broad (FWHM ~ 2850 km s" 1 ) components in the spectrum of RX J1334.2+3759. 
The narrow-line region of RX J1334. 2+3759 appears to be significantly different from 
those of normal Seyfert galaxies. A possible explanation for the observed properties 
of the narrow line region and the broad line region is suggested in terms of density 
, enhancements. 
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1 INTRODUCTION 

Narrow-line Seyfert 1 (NLS1) galaxies are considered to be 
a special class of "normal" Seyfert 1 galaxies because of 
their peculiar properties that distinguish them from the lat- 
ter class. They are characterized by their optical spectra 
having permitted lines that are narrower than those in the 
normal Seyfert 1 galaxies, e.g., full width at half maximum 
(FWHM) of H/3 line is <5 2000 km s"\ relatively weak for- 
bidden lines, l°Iin**°°l < 3 (Osterbrock & Pogge 1985; 
Goodrich 1989), and strong Fe II emission. X-ray observa- 
tions with Einstein (e.g., Puchnarewicz et al. 1992) and 
ROSAT (e.g., Boiler, Brandt, & Fink 1996) have revealed 
that NLS1 galaxies have distinctive soft X-ray properties as 
well. NLS1 galaxies show steep soft X-ray spectra with little 
or no absorption above the Galactic values along their re- 



spective lines of sight (Grupe et al. 1998). They often show 
rapid and large amplitude as well as long-term X-ray vari- 
ability (Boiler et al. 1993, Brandt, Pounds, & Fink 1995, 
Grupe et al. 1995a, b). In spite of the dominance of soft X- 
ray emission, the soft X-ray luminosities of NLS1 galaxies 
are similar to those of normal Seyfert Is. High luminosity 
analogs of the NLS1 class, however, do exist. Proto- types of 
this class being - a narrow-line QSO (or NLQSO) I Zw 1 
(Phillips 1976), PHL 1092 (Brandt 1995; Forster & Halpern 
1996; Lawrence et al. 1997), PKS 0558-504 (Remillard et 
al. 1991), and PG 1404+226 (Ulrich et al. 1999). 

The spectral energy distribution (SED) from far- 
infrared (FIR) to X-rays of NLS1 galaxies appears to be 
similar to that of broad-line Seyfert 1 galaxies, but the UV 
luminosities of NLS1 galaxies tend to be smaller than those 
of Seyfert Is (Rodriguez-Pascual, Mas-Hesse, & Santos- 
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Lleo 1997). The lower UV luminosity of NLS1 galaxies 
compared to normal Seyfert Is could be due the shift of 
the Big Blue Bump (BBB) towards higher energies. The 
steep soft X-ray spectrum could be the high energy tail of 
the BBB (Mathur 2000). To explain the relatively narrow 
width (FWHM~ 2000 km s" 1 ) of the H/3 line in the NLS1 
galaxies as compared to those (FWHM~ 5000 km s _1 ) in 
the Seyfert 1 galaxies, Wandel (1997) has argued that the 
steeper soft X-ray continuum has a higher ionizing power, 
resulting in an extended broad line region (BLR) , and hence 
a smaller velocity dispersion and narrower emission lines. Al- 
ternative views involve the virialized motion of similar size 
BLR around a smaller black hole mass in NLSls compared 
to those in Seyfert Is (see Laor et al. 1997a; Brandt & Boiler 
1998), and disc inclination effects (e. g Osterbrock & Poggc 
1985; Goodrich 1985). The weakness of the forbidden line, 
[O III]A5007, relative to H/3 has not been understood. The 
importance of the decomposition of the Balmer lines, Ha 
and H/3, into narrow and broad components has been real- 
ized only recently in order to study the narrow line region 
(NLR) and BLR emission line ratios in the NLS1 galax- 
ies. Recent work of Rodriguez- Ardilla et al. (2000) has re- 
vealed that emission line ratios from the Narrow Line Region 
(NLR) of NLS1 are different from those observed in Seyfert 
galaxies, e.g. [ °"^ 5U " 7 = 0.8 - 5.0 in the NLR of 7 NLS1 
galaxies, which is significantly different from the value ~ 10 
observed in the NLR of Seyfert galaxies. The reason for the 
observed difference has not been understood satisfactorily. 
Also, there are only a few NLS1 galaxies for which the NLR 
omission line ratios have been determined. Therefore, it is 
important to study the NLR and BLR emission line ratios 
in the other NLS1 galaxies. 

In this paper, we present soft X-ray and optical emis- 
sion line properties of a new narrow-line QSO (NLQSO) - 
RX J1334.2+3759, and investigate the origin of the NLR and 
BLR emission line ratios. The object, RX J1334. 2+3759, 
was identified as an ultra-soft X-ray source by Singh et 
al. (1995) based on WGACAT (White, Giommi, & Angelini 
1994). Using deep ROSAT PSPC observations, the opti- 
cal counterpart of the X-ray source RX J1334. 2+3759 was 
identified by McHardy et al. (1998). The basic parameters of 
RX J1334. 2+3759 are given in Table 1. The paper is struc- 
tured as follows. In the next section, we describe the X-ray 
and the optical spectroscopic observations followed by anal- 
yses of the X-ray data and the optical spectrum in §3. We 
discuss our results in §4 followed by conclusions in §5. 

Throughout the paper, luminosities are calculated as- 
suming isotropic emission, a Hubble constant of Ho = 
75 km s _1 Mpc -1 and a deceleration parameter of qo — 
unless otherwise specified. 

2 OBSERVATIONS 
2.1 X-ray 

The 

region of the sky containing the source RX J1334. 2+3759 
was observed twice with the ROSAT (Triimper 1983) Posi- 
tion Sensitive Proportional Counter (PSPC) during 1991- 
1993, and once with the High Resolution Imager (HRI) 
(Pfeffermann et al. 1987) in 1997 June 19-July 16. The 



two PSPC observations carried out during 1991 June 23- 
26, and 1993 June-July together comprise the second deep- 
est ROSAT PSPC survey. These observations were tar- 
geted at a(2000) = 13 h 34 m 37.0 s , 5(2000) = +37° 54' 44" 
in the sky with an extremely low obscuration due to mat- 
ter in our Galaxy (N H = 7.9 x 10 19 cm" 2 ). The PSPC 
deep survey data has already been reported by McHardy 
ct al. (1998). The HRI observation was also a deep survey 
(exposure time = 201513 s) and was carried out in the same 
region of sky as the PSPC deep surveys. The details of the 
ROSAT observations are given in Table 2. The off-sets of 
RX J1334. 2+3759 from the field centres are also listed in 
Table 2. ROSAT X-ray data corresponding to the above 
observations were obtained from the public archives main- 
tained at the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC) in USA, and analyzed by us us- 
ing PROS and FTOOLS. 



2.2 Optical 

An optical spectrum of RX J1334. 2+3759 was obtained on 
the night of 1994 April 7 with the Multi Object Spectro- 
graph (MOS) on the 3.6-m Canada France Hawaii Telescope 
(CFHT) as a part of a program to identify X-ray sources 
(McHardy et al. 1998). A 300 1 mm" 1 grism in the first or- 
der with a Lorel3 CCD detector was used to cover a wave- 
length range of 4000 A-9000 A with ~ 15 A resolution. The 
integration time was 2700 s. For details of the reduction see 
McHardy et al (1998). The flux calibration is somewhat un- 
certain because the slitlets were not aligned at the parallac- 
tic angle. 

Optical R band image of RX J1334. 2+3759 was ob- 
tained from the 2.5-m Isaac Newton Telescope (INT) using 
the wide-field camera on the night of 1999 April 19. The ex- 
posure time was 1200 s. The pixel scale is 0.333 arcsec per 
pixel. The seeing (FWHM) was ~ 1.3". 



3 ANALYSIS & RESULTS 

The X-ray source, presented in this paper, was identified 
with an optical object (R= 19.5) by overlaying the contours 
of high resolution (central full width half maximum of ~ 4") 
X-ray image obtained from ROSAT HRI observations onto 
the R-band optical image. HRI images were extracted and 
smoothed by convolving with a Gaussian of a = 2" using the 
PROS software package. X-ray contours overlaid onto the R- 
band optical image are shown in Figure 1. Two fainter point 
objects about 4 and 10 arcsec away from RX J1334. 2+3759 
are seen on the overlay. The R magnitudes of these objects 
are ~ 22.3 (nearest; hereafter object A) and ~ 22.9 (object 
B) estimated from the R band image and known R magni- 
tude of RX J1334.2+3759. Object B is unlikely to be the 
optical counterpart of the X-ray source because of its dis- 
tance from the X-ray peak. On the other hand, object A is 
unresolved and could contribute the X-ray intensity. In order 
to investigate the amount of X-ray emission from the object 
A, we calculate the ratio of soft X-ray to R-band fluxes, -y^, 
for 22 QSOs including RX J1334. 2+3759, which are within 
a circle of radius 15 arcmin (Dewangan et al. 2001). The flux 
ratios were calculated using the relation 
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Table 1. Basic parameters of RX J1334.2+3759 



Position 1 : a(J2000) = 13 h 34 m 10.6 s ; <5(J2000) = +37° 59' 56.0". 
Redshift 2 : z = 0.3858 ± 0.0004. 
Magnitude 1 : R = 19.5. 

iMcHardy ct al. (1998) 
2 Present paper 



Table 2. Details of ROSAT observations of RX J1334.2+3759 



Serial 


Sequence 


Instrument 


Offset 


Start Time 


End Time 


Exposure 


Count Rate" 


No. 


No. 




arcmin 


Y, M, D, UT 


Y, M, D, UT 


Time (s) 


10~ 2 cnt s- 1 


1. 


RH900717N00 


HRI 


8.207 


1997 06 04 16:12:58 


1997 07 13 22:26:43 


201513 


0.41 ± 0.02 


2. 


RP900626N00 


PSPC 


8.207 


1993 06 19 22:24:46 


1993 07 16 23:06:28 


37658 


1.93 ± 0.07 


3. 


RP700283N00 


PSPC 


8.207 


1991 06 23 20:59:40 


1991 06 26 20:44:16 


71803 


1.71 ± 0.09 



a Mean count rates after background subtraction in the energy band of 0.1-2.4 keV. 



log(^)=log(f x ) + r -^+3.5 (1) 

JR 2.5 

where m_R is the R-band magnitude. The soft X-ray fluxes 
of all the QSOs have been derived from the ROSAT PSPC 
observation of 1991 listed in Table 2 (Dewangan et al. 2001). 
The ratio, s -f-, for the QSOs ranges from 0.003 to 0.040 with 
a mean value of 0.014 ± 0.010. If we assume the object A 
to be the counterpart of the X-ray source, then the ratio 
(y^-) is estimated to be ~ 0.20 for A. A comparison of this 

ratio with that obtained for the QSOs shows that j^- for the 
object A is much higher than that for the QSOs. However, 
^ for RX J1334.2+3759 is 0.017 which is similar to that 
of other QSOs. If the object A is also similar to the QSOs, 
the contribution of this source to the total X-ray emission 
is expected to be ~ 7%. 

The total source counts for RX J1334. 2+3759 were ob- 
tained from the unsmoothed PSPC images using a circle 
of radius of 1.6' centered on the peak position, and af- 
ter subtracting the background estimated from five nearby 
circular regions with their centres ~ 7.0' away from the 
source. This was done using the xselect program in the 
FTOOLS (version 5.0) software package. The HRI count 
rate for RX J1334. 2+3759 was extracted from a circle of 
radius 25" centered on the peak position, after subtract- 
ing the background estimated from an annulus of width 1.5' 
with the same centre and with an inner circle radius of 2'. 
The count rates thus estimated are given in Table 2. 



3.1 X-ray Light Curves 

In order to investigate the time variability of soft X-ray 
emission from RX J1334. 2+3759, we have extracted the 
light curves from the ROSAT PSPC observations. The light 
curves for the source and the background were extracted us- 
ing the 'xselect' package in the PSPC energy band of 0.1-2.4 
keV containing all the X-ray photons falling within "good 
time intervals". The time bin sizes are 500 s so that on an 
average each bin has ~ 10 counts. The source regions and 
the background regions were the same as described above. It 
was found that the background was highly variable during 
the observation of 1993 June 19. Therefore, the light curve 



of RX J1334. 2+3759 obtained from the 1993 observation is 
not suitable for variability studies. During the observation 
of 1991 June 23, the background was reasonably constant. A 
constant count rate fit to the background light curve gives 
the best-fit minimum value of \ 2 ~ 73.54 for 75 degrees of 
freedom. The background subtractions were carried out af- 
ter appropriately scaling the background light curve to have 
the same area as the source extraction area. The background 
subtracted light curve of RX J1334. 2+3759 and the back- 
ground light curve are shown in Figure 2. A remarkable vari- 
ability in the soft X-ray flux from RX J1334. 2+3759 can be 
seen in Fig. 2. A constant count rate fit to the light curve of 
RX J1334. 2+3759 gives minimum value of \ 2 of 154.6 for 75 
degrees of freedom. X-ray emission from RX J1334. 2+3759 
changed by a factor of ~ 2 on time scales of 20000 — 40000 s. 
Variability on shorter time scales is also observed on several 
occasions notably at the beginning of the observation, af- 
ter ~ 200000 s, and after ~ 240000 s. The most significant 
and extreme variable event is that observed around 200000 s 
which is shown in Figure 3 with a smaller bin size of 200 s. 
The 9 data points to the left of the variability event (only 
two data points are shown in Fig. 3) have a mean value of 
(2.21 ± 1.26) x 10~ 2 count s" 1 , we treat this average count 
rate as the rate in the quiescent state from which the flare 
arises. After 200200 s from the beginning of the observation, 
the PSPC count rate began to increase reaching a maximum 
of (8.90 ± 2.25) x 10~ 2 count s _1 at 200600 s. Thus, an in- 
crease in the count rate by a factor of ~ 4 in < 400 ± 141 s 
(< 380 ± 134 s in the rest frame) is detected. Thus, for 
the variability event shown in Fig. 3, we find a change 
in the PSPC count rate of (6.69 ± 2.58) xl0~ 2 count s" 1 . 
This corresponds a change in the intrinsic luminosity of 
(7.43 + 2.86) x 10 44 erg s _1 in the energy band of 0.1-2.0 keV 
in the rest frame time interval of < 380 ± 134 s. 

In order to investigate the long term (time scale of a few 
years) variation in the X-ray intensity of RX J1334. 2+3759, 
we have converted the HRI count rate, obtained from the 
observation of 1997 June 4, into equivalent PSPC count rate 
using the best-fit model parameters (Nh = 3.3 x 10 20 cm~ 2 , 
Fx = 3.8) obtained from the joint-fit to the PSPC spectra 
(see §3.2). The equivalent PSPC count rate is estimated to 
be (2.22 ± 0.11) x 10~ 2 count s~ x , which is similar to the 
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Figure 1. Contours of ROSAT HRI intensity of RX J1334. 2+3759 overlaid on the R-band optical image. For clarity, only the contours 
at 6%, 10%, 20%, 30%, 40%, 50%, 60% of the peak intensity are shown. The X-ray contours have been generated from the HRI image 
after smoothing by a Gaussian of <r = 2". The HRI image was created from the data observed on 1997 July 7 (exposure time = 201513 s). 



PSPC count rates obtained during the observations of 1991 
and 1993. 

3.2 X-ray Spectral Analysis 

Photon energy spectra of RX J1334. 2+3759 were accu- 
mulated from their PSPC observations shown in Table 2. 
The same regions for the source and the background, as 
stated above, were used. The ROSAT PSPC pulse height 
data obtained in 256 pulse height channels were appropri- 
ately re-grouped to improve the statistics. X-ray spectra of 
RX J1334. 2+3759 from the two observations thus obtained 
are shown in Figure 4. 

We used the XSPEC (Version 11.0) spectral analysis 
package to fit the data with spectral models. This requires 
a knowledge of the response of the telescope and the de- 
tector. An appropriate response matrix, provided by the 
ROSAT GOF at HEASARC, was used to define the energy 
response of the PSPC. The off-axis calibration of the tele- 
scope has been appropriately taken into account by creating 
and using an auxiliary response. The ROSAT PSPC spec- 
tra of RX J1334. 2+3759, shown in Figure 4, were used for 



fitting spectral models. The spectra from the two observa- 
tions were first fitted separately with redshifted power-law 
model with photon index, Tx, and absorption due to an 
intervening medium with the absorption cross-sections as 
given by Balucinska-Church & McCammon (1992) and us- 
ing the method of x 2 -minimization. The results of these fit- 
ting and the best-fit model parameters are shown in Table 
3. The errors quoted, here and below, were calculated at the 
90% confidence level based on Xmm+2-71. As can be seen 
in Table 3, the best-fit model parameters derived from the 
two observations are similar within errors. In order to bet- 
ter constrain the model parameters, we have jointly fitted 
spectral models to the two observed spectra and below, we 
discuss the joint model fitting in detail. The results of this 
fitting and the best-fit spectral model parameters are listed 
in Table 3. The simple model, (power-law + absorption at 
z = 0), with Fx in the range of 3.4-4.0 and Nh in the range 
of 1.5 x 10 20 - 2.5 x 10 20 cm" 2 is a good fit to both the spec- 
tra, as evidenced by the minimum reduced x 2 (xt) value of 
0.95 for 43 degrees of freedom. The absorbing column den- 
sity thus derived by fitting the power-law model is in excess 
of the Galactic value (Nh = 7.9 x 10 19 cm -2 ) measured 
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Figure 2. Background subtracted ROSAT PSPC light curve of RX J1334.2+3759 observed on 1991 June 23 and the simultaneously 
obtained background light curve (solid line). The bin sizes in the light curves is 500 s. The start time of the observation was 1991 June 
23 20:59:40 UT. 



from 21-cm radio observations along the line of sight to the 
source (Dickey & Lockman 1990). This indicates that all the 
X-ray absorption may not be only due to matter in our own 
Galaxy but also due to matter local to the source. We next 
fixed the absorbing column to the Galactic value and fitted 
the power-law model to both spectra. There is an increase 
in the xi value and the fit becomes worse (xt = 1-57 for 44 
degrees of freedom) indicating that an additional component 
is needed. In order to determine the absorbing column local 
to the source, we introduced an additional absorbing col- 
umn at the source redshift (z = 0.3858) and carried out the 
power-law fitting. The best-fit spectral parameters are again 
given in Table 3. The intrinsic absorption column density 
(AN H ) thus derived is indeed small (3.3tH x 10 20 cm" 2 ). 
The best- fit photon index is 3.81q'|. 

Based on the xt values, it is clear that a redshifted 
power-law model with Galactic plus intrinsic Nh is pre- 
ferred over the model with ANh = 0. The power-law 
model absorbed by the Galactic and the intrinsic column 
is a significant improvement over the power-law model ab- 
sorbed by the Galactic column alone, at 99.99% confi- 
dence level based on an F-test (F-statistics value = 29.7 
and probability = 2.28 x 10" 6 ). Allowed ranges of Tx 
and ANh are shown in Figure 5. The observed X-ray 
flux from the source RX J1334. 2+3759 is estimated to be 
8.7xl0~ 14 erg cm -2 s _1 in the energy band of 0.1-2.0 keV. 



The intrinsic flux , estimated by setting Nh to zero, is found 
to be 5.4 x 10~ 13 erg cm~ 2 s _1 in the energy band of 0.1-2.0 
keV. The intrinsic X-ray luminosity in the energy band of 
0.1-2.0 keV is calculated to be 2.2 x 10 44 erg s" 1 . We have 
also fitted redshifted blackbody models absorbed by an in- 
tervening medium. The absorbing column density was fixed 
to the Galactic value as otherwise unphysically low values of 
Nh , lower than the Galactic Nh , are obtained. The best-fit 
model parameters are given in Table 3. It is clear from the 
xl that the PSPC spectra are not well represented by a sim- 
ple blackbody model, however, a temperature of ~ 135 eV 
is inferred for the blackbody. 

3.3 Optical Spectrum 

The optical spectrum of RX J1334. 2+3759 is shown in 
Figure 6. The signal-to-noise ratio of the spectrum is ~ 
13, as measured from the dispersion in the continuum 
region 6050 A — 6150 A in the spectrum. Strong emis- 
sion lines of Balmer Ha, H/3, and the forbidden line [O 
III]A5007 are readily observed. To measure the redshift of 
RX J1334. 2+3759, we fitted Gaussians to the upper half of 
the Ha, H/3, and [O III]A5007 line profiles. Using the fitted 
positions and the rest wavelengths of these lines, the red- 
shift of RX J1334.2+3759 is derived to be 0.3858 ± 0.0004. 
The redshift quoted here is the average value derived from 
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Figure 3. The most extreme variability event observed from RX J1334. 2+3759 about 2 X 10 5 s after the beginning of the observation 
of 1991 June 23. The bin sizes in the light curve is 200 s. The light curve is background subtracted. The start time of the observation 
was 1991 June 23 20:59:40 UT. 



the above three lines and the error is the dispersion. The 
spectrum, shown in Fig. 6, is corrected to the rest wave- 
lengths. Apart from the emission lines mentioned above, we 
have also identified several other emission lines such as [Ne 
V]A3426, [O II]A3727, [Ne IHJA3869, He A3969, H 7 A4340, 
H/3, Fe IIA4924, and He IA5876. Also present is the Fe II 
emission seen as two broad humps at 4450 A-4700 A and 
5150 A-5350 A in the spectrum. All the above lines have 
been marked in Fig. 6. 

The Fe II blends complicate the measurement of line 
strengths of H/3, [O III]A5007 etc. due to their contamina- 
tion. In order to reliably estimate their strength and to mea- 
sure the line fluxes, we have adopted the method of Boro- 
son & Green (1992) and used their Fe II template. In order 
to use the Fe II template, first we ensured that the tem- 
plate and the spectrum of RX J1334. 2+3759 are correctly 
redshifted and refer to the rest wavelengths. The template 
was then broadened to the FWHM of H/3 line in the spec- 
trum of RX J1334. 2+3759 by convolving with a Gaussian 
of (j = 7.5 A. The smoothed Fe II template was scaled ap- 
propriately to match the intensity of the observed Fe II and 
then subtracted from the spectrum of RX J1334. 2+3759. 
The procedure of Fe II subtraction is depicted in Figure 7. 
Note that the Fe II subtracted spectrum has been shifted 
downwards for clarity. As a result of the Fe II subtraction, 
the continuum around the H/3 line has flattened implying 
that the Fe II model fits well the observed Fe II in the spec- 
trum of RX J1334. 2+3759. The rest frame equivalent width 



of Fe II was determined from the flux in the scaled Fe II 
template between the rest wavelengths 4250 A and 5880 A 
and the continuum flux density at 5050 A in the Fe II sub- 
tracted spectrum. The equivalent width of Fe II is listed in 
Table 4. 

In order to measure the parameters of strong lines other 
than Fe II, we have used the continuum subtracted and 
Fe II corrected spectrum of RX J1334.2+3759. The con- 
tinuum from the Fe II subtracted spectrum was removed 
by fitting a low-order polynomial after excluding the emis- 
sion lines. We have assumed that the emission lines can be 
represented by a single or a combination of Gaussian pro- 
files. As a first step, we fitted a single Gaussian profile to 
the H/3 line. The FWHM velocity of the best-fit Gaussian is 
1557 km s" 1 . It was, however, found that the broad wings as 
well as the peak of the H/3 line are not well fitted by a single 
Gaussian profile. It was necessary, therefore, to introduce 
an additional component in order to represent the observed 
profile adequately. We then fitted two Gaussian profiles si- 
multaneously to the H/3 line. The decomposition of the H/3 
line into two Gaussians is shown in Figure 8. The positions 
of the Gaussians were fixed at the rest wavelength of H/3 
line while the amplitudes and the FWHMs were varied to 
achieve the best fit. The residuals from the best fitted pro- 
file were found to be quite similar to the noise level near 
H/3. Thus we find that two Gaussians, a narrow component 
(FWHM ~ 877 km s" 1 ) and a broad component (FWHM 
~ 2853 km s _1 ) describe the observed H/3 profile much bet- 
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Figure 4. Spectral data from 2 observations and a model (histogram) fitted simultaneously to the ROSAT PSPC X-ray spectra of 
RX J1334. 2+3759. The fitted model is a rcdshifted power-law modified by absorbing columns present in our own Galaxy as well as in 
RX J1334.2+3759. 



ter than a single Gaussian profile. The parameters of the 
H/3 line obtained from the best fitting profile are listed in 
Table 4. The equivalent widths, given in Table 4, refer to the 
continuum flux density, obtained from continuum fitting, at 
the rest wavelength of H/3. The [O IIIJA5007 line has not 
been resolved. Its width is found to be comparable to the 
instrumental resolution. The strength of [O III]A5007 line 
was determined by integrating the line profile. 

Due to the poor instrumental resolution (FWHM 
~ 15 A), it was not possible to decompose Ha and 
[N II]AA6548,6583 lines. To measure reliably the FWHM 
and the equivalent width of Ha, it is necessary to correct 
for the contribution of [N II] AA6548, 6583. For this pur- 
pose, we have created templates for the [N II]AA6548, 6583 
lines from the [O III]A5007 line. These lines are expected 
to have similar widths and are unresolved in the spectrum 
of RX J1334.2+3759. First, we isolated the [O III]A5007 
line in the wavelength region 4992 A-5018 A from the con- 
tinuum subtracted and the Fe II corrected spectrum of 
RX J1334.2+3759. The [O IIIJA5007 line was then shifted to 
the rest wavelength of [N II]A6583. To get the template for 
[N IIJA6583, the shifted [O IIIJA5007 line was scaled to have 
flux that was 35% of the actual flux of the [O III]A5007 line, 
since the contribution of [N II]A6583 to Ha is ~ 35% of the 
flux of [O III]A5007 (Ferland & Osterbrock 1986). We did 
not change the FWHM of [N II]A6583 template from that 



of [O III]A5007. To make the template for the [N II]A6548 
line, the [N IIJA6583 template was shifted to the rest wave- 
length of [N IIJA6548 and then scaled so that the ratio of 
fluxes of [N IIJA6583 and [N II]A6548 is the same as the 
theoretical ratio of 2.96. To correct for the contribution of 
[N II]AA6548, 6583 to the Ha line, the two templates, created 
above, were subtracted from the continuum subtracted and 
Fe II corrected spectrum of RX J1334.2+3759. This subtrac- 
tion has little effect on the Ha profile because the intensities 
of [N II] AA6548, 6583 lines are much smaller than that of Ha. 
This is expected since the flux ratio, ^° II ^ o0m is very small 
(~ 0.14) for RX J1334.2+3759 which gives the flux ratio, 
[at /i"]AA6o48,6j83 ^ § q\q assuming the case B value of the 
flux ratio ~ 3.1 (e.g. Osterbrock 1989). The Ha profile is 
again poorly fitted by a single Gaussian. Applying the same 
decomposition technique, as was done for the H/3 profile, we 
find that the profile of Ha is best fitted by two Gaussians, 
one with a narrow component (FWHM ~ 883 km s _1 ) and 
another with a broad component (FWHM ~ 2843 km s _1 ). 
The two component profile fit to the Ha is shown in Figure 
9. The equivalent widths, listed in Table 4, refer to the flux 
density at the rest wavelength of Ha in the Fe II corrected 
spectrum. 
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Figure 5. Allowed ranges of power-law photon index and 'excess' Njj (ANh ) (at z=0.38) for 68%, 90%, and 98% confidence based on 
counting statistics derived from the joint-fit of the 2 PSPC spectra for RX J1334. 2+3759. The '+' marks the best-fit value. 



4 DISCUSSION 

RX J1334. 2+3759 is luminous in optical and soft X-rays. 
Its optical luminosity (Mr ~ —22.7 assuming Ho = 
50 km s _1 Mpc -1 ) and soft X-ray luminosity (Lx — 2.8 x 
10 44 erg s _1 in the energy band of 0.1 — 2.0 keV) are sim- 
ilar to that of a low luminosity QSO. The object has been 
classified as a QSO by McHardy et al. (1998). Based on op- 
tical spectroscopy and analysis of publically available X-ray 
data from ROSAT observations, we find that the newly dis- 
covered QSO RX J1334.2+3759 belongs to the NLS1 class. 
Therefore, it is referred to as NLQSO here. Its X-ray and 
optical characteristics are further discussed below. 



4.1 Soft X-ray Variability 

Soft X-ray emission from RX J1334. 9+3759 is highly vari- 
able (see Fig. 2 & 3). During the 1991 observations, 
RX J 1334.9+3759 showed soft X-ray variability on times 
scales of ~ 20000-40000 s by a factor of ~ 2. Rapid variabil- 
ity events have also been detected from RX J1334. 9+3759. 
The most significant and extreme variable event, shown in 
Fig. 3, has > (1.95 ± 1.02) x 10 42 erg s~ 2 , which is sim- 
ilar to that of the extreme variable event observed from 
PHL 1098 (Brandt et al. 1999). Straightforward application 
of the efficiency (rj) limit: ij>4.8x 10~ 43 §| (Fabian 1979), 
results in an extremely high efficiency (r) > 0.93 + 0.49). This 
can be compared with r) ~ 0.3 for an optimally accreting 
Kerr BH rotating at the maximum plausible rate (Thornc 
1974). Thus the radiative efficiency for the extreme variable 
event shown in Fig. 3 appears to be substantially larger and 
suggests relativistic effects responsible for the extreme vari- 



able event. The suggestion that relativistic effects may be 
responsible for the extreme variability event can also be in- 
ferred as follows. The variability time scale, in the absence 
of relativistic effects, provides an upper limit to the size of 
the X-ray emitting region, R < cAt. For RX J1334.2+3759, 
R < 1.14 x 10 13 cm. SEDs of a sample of NLS1 galaxies 
(Rodriguez-Pascual, Mas-Hesse, & Santos-Lleo 1997) show 
that the soft X-ray luminosity of NLSls can, at most, be 9% 
of the bolometric luminosity. If we assume that the SED of 
RX J1334.2+3759 is similar to that of the other NLSls, then 
the bolometric luminosity of RX J1334. 2+3759 is expected 
to be <; 2.4 x 10 45 erg s _1 . If this luminosity is within the 
Eddington limit, then RX J1334. 2+3759 must have a cen- 
tral mass, M > 2 x 10 7 M0. In that case, the radius of the 
last stable orbit, 3Rs > 1.8 x 10 13 cm, would then be larger 
than the size of the X-ray emitting region, so photons should 
not escape. 

Extremely rapid X-ray variability has also been ob- 
served in three other NLSl-class objects: PKS 0558-504 
(77 > 1.5; Remillardet al. 1991), IRAS 13224-3809 (rj > 0.09; 
Boiler et al. 1997), and PHL 1092 (r) > 0.6; Brandt et 
al. 1999) and relativistic effects have been suggested to be 
responsible for the enhanced X-ray variability. 

4.2 Soft X-ray Spectral characteristics 

The soft X-ray spectrum of RX J1334. 2+3759 is very 
steep and is well represented by a power-law of photon 
index, Fx = 3.8^q 3 with an excess absorption, ANh ~ 

3.3 x 10 20 cm~ 2 , local to the source and apart from the 
Galactic absorption, N H ~ 7.9 x 10 19 cm~ 2 (see Table 
3). Thus the soft X-ray spectrum of RX J1334.2+3759 is 
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Table 3. Best-fit Model spectral parameters for RX J1334.2+3759 



Data 


Model" 


N H 

(10 19 cm- 2 ) 


A N H b 
(10 20 cm- 2 ) 


Fx 

or kT (eV) 


fx 


L x d 




1991 (June 23-July 26) 


pi + abs Gal 


18.112;? 




3 6+ ' 3 


8.8 


1.71 


0.97/24 




pi + abs Gal 


7.9 (fixed) 




o 1+0.1 

^•-'■-o.i 


8.22 


0.72 


1.58/25 




pi + abSGal + abs S ource 


t .y ^nxcQ ) 


9 8+ 1 - 7 


n 7 +0.4 


o.o 


1 Q 
i .y 


u.yo / 


1993 (June 19- July 16) 


pi + abs Gal 






4 o+ ' 8 


8.55 


2.94 


0.960/16 




pi + abs Gal 


7.9 (fixed) 




2 9+ ' 2 


8.3 


0.68 


1.64/17 




pi + abs Gal + abs a0 u rC e 


7.9 (fixed) 




4 2+ 10 


8.6 


3.66 


0.96/16 


1991 and 1993-f 


pi + abs Gal 


19.8+11 




o 7 +0.3 
°- ' -0.3 


8.7 


1.90 


0.95/43 




pi + abs Ga i 


7.9 (fixed) 




3 n+ ' 1 


8.20 


0.33 


1.57/44 




pi + abs Gal + abs sourco 


7.9 (fixed) 


Q Q + 1.6 


o O+0.3 
J - S -0.3 


8.7 


2.18 


0.95/43 




abs + blackbody 


7.9 (fixed) 




135 


9.73 


0.4 


1.74/44 



a pl is rcdshiftcd simple power-law model. abs Ga i(abs SOU rcc) is the photo-electric absorption model at z = 0(z = 0.3858) using 
Balucinska-Church, and McCammon (1992) cross-sections, blackbody is the redshiftcd blackbody model. 
'Excess Nh over the Galactic value. 

c Observed flux in units of 10- 14 erg cm -2 s _1 in the energy band of 0.1 — 2.0 keV. 
^Intrinsic soft X-ray luminosity in units of 10 44 erg s _1 in the energy band of 0.1 — 2.0 keV. 
e Minimum reduced \ 2 for v degrees of freedom. 
/Data are jointly fit to the same spectral models. 



Table 4. Optical emission line parameters of RX J1334. 2+3759 



Emission line 


Broad Component 


Narrow Component 




FWHM 


Equivalent Width 


FWHM 


Equivalent Width 




km s~ 1 


A 


km s _1 


A 


H/3 


2853 ± 520 


15.8 ±3.6 


877 ± 91 


19.5 ±3.8 


[O III]A5007 








5.0 ±0.5 


Ha 


2843 ± 479 


33.8 ±7.5 


883 ± 41 


54.5 


Fe II 




159 







steeper than those of normal Seyfert Is [< Yx > (90% 
range)= 2.0 — 2.7], and similar to those of NLS1 galaxies 
[Tx(90% range) = 2.3 - 3.7] (Grupe et al. 1998). The in- 
trinsic absorption, inferred from the power-law model fit to 
the spectra of RX J1334.2±3759 is not high, and is simi- 
lar to those found in normal Seyfert Is and NLS1 galaxies. 
The steeper power-law index and blackbody model fit to 
the PSPC spectra of RX J1334.2±3759 indicate the ultra- 
soft nature of this object. The derived temperature of the 
blackbody, kT ~ 135 eV, is similar to those found in NLS1 
galaxies (Brandt & Boiler 1998). That this object is simi- 
lar to the NLS1 galaxies in its nature, is further supported 
by the optical spectroscopic properties of RX J1334. 2+3759 
(see below). The excess soft X-ray emission of NLS1 galaxies 
is usually attributed to a higher accretion rate compared to 
the Eddington accretion rate (m = ^ ) (Pounds, Done, 



& Osborne 1995; Brandt & Boiler 1998). The bolometric lu- 
minosity of RX J1334.2+3759, L bol > 2.4 x 10 45 erg s" 1 , is 
about a factor of ~ 2 higher than the Eddington luminosity 
for a 10 7 M BH. Thus Eddington or super-Eddington ac- 
cretion rate is required in RX J1334. 2+3759 provided that 
the mass of the central BH is ~ 10 7 Mq. 



4.3 Optical Spectral Characteristics 

The optical spectrum of RX J1334. 2+3759, shown in Fig. 6, 
appears to be typical of NLS1 galaxies. The FWHM veloc- 
ity (~ 2850 km s _1 ) of the broad components of Balmer 
lines in the spectrum of RX J1334. 2+3759 is similar to 
that found in the other NLS1 galaxies but is significantly 
lower than in the Seyfert 1 galaxies. For example, Grupe et 
al. (1999) found mean values of the broad components of 
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Figure 6. Optical spectrum of RX J1334. 2+3759. Vertical scale is the relative flux. The spectrum has been smoothed by a box filter of 
width 12.5 A. The absorption feature at 5500 A is due to atmospheric absorption. 



H/3 to be 2790 ± 160 km s" 1 for a sample of NLSls, and 
4210 ± 360 km s _1 for a sample of Seyfert Is. Similarly, 
the ratio of fluxes of narrow and broad components of Ha 
and H/3, f|a = i.g ± 0.46, and = 1.26 ± 0.37, for 

RX J1334.2+3759 are similar to those 'found in other NLS1 
galaxies (Rodriguez- Ardila et al. 2000). This implies that 
the relative contribution of the broad components to the 
line flux is greatly reduced in RX J1334. 2+3759 compared 
to that in Seyfert 1 galaxies for which the ratio of narrow to 
broad component is around 0.1. 

The narrow components of the permitted lines and 
the forbidden lines such as [O III]A5007, [N II]A6583 
etc. arc thought to arise from the NLR of active galac- 
tic nuclei (AGNs). For RX J1334. 2+3759, the flux ratio 
of [O III]A5007 line and the narrow component of H/3, 

^° I h0 X6007 i ls O'^O ^ 0-06 which is much smaller than the 
ratio of ~ 10 found in Seyfert 2 galaxies. NLS1 galaxies, on 
the other hand, show values ranging from 0.5 to 5 for the 
l ° "h'IT' 07 ratio ( see Rodriguez- Ardila et al. 2000). Thus, 
NLR of RX J1334+3759 is somewhat different from the 7 
NLS1 galaxies in the sample of Rodriguez- Ardila et al. 2000. 
RX J1334.2+3759 shows strong Fe II emission. The 



ratio of fluxes of Fe II and H/3, 



H0 



4.5 for 



RX J1334. 2+3759. This ratio is similar to those found in 
the other NLS1 galaxies (see Grupe et al. 1999). 

In the following, we try to explain the observed flux ra- 
tios from the NLR and BLR of RX J1334. 2+3759 and the 



other NLS1 galaxies in terms of density enhancements. In 
the same picture, we also try to explain the well known anti- 
correlation between the slope of the soft X-ray continuum 
and the FWHM of H/3 line found in NLS1 galaxies. Ed- 
dington or Super-Eddington accretion rates, thought to be 
responsible for the steep soft X-ray continua in NLSls (see 
Pounds, Done, & Osborne 1995; Brandt & Boiler 1998), may 
result in outflows. In the NLR of NLS1 galaxies, the flux ra- 



tio 



[O II I] A5007 
H/3„ 



is smaller than the ratio found in the normal 



Seyfert galaxies. This smaller value can be produced by a 
smaller size of the NLR in NLS1 galaxies. Outflows can en- 
hance the density of the BLR and extend its size outwards, 
thereby reducing the size of the NLR. In this scenario, the 
flux of the [O III]A5007 line will be reduced while that of the 
narrow component of permitted lines will not be. Rather, 
the flux in the narrow components of the permitted lines 
is expected to increase due to an enhancement of the den- 
sity. Thus, a smaller ratio of [O III]A5007 and H/3 n can be 
produced in the NLR of NLS1 galaxies. 

The FWHM of the broad component of H/3 line is sig- 
nificantly smaller in the NLS1 galaxies than those found in 
Seyfert 1 galaxies. It seems that in the NLS1 galaxies, ap- 
propriate physical conditions, where the maximum number 
of broad component H/3 photons can be produced, are lo- 
cated at a larger radial distance from the nucleus than in 
the normal Seyfert 1 galaxies. This radial shift seems to 
be possible due to outflows and a resulting enhancement 



© 0000 RAS, MNRAS 000, 000-000 



Optical and X-ray characteristics of a newly discovered narrow-line QSO : RX 3133^.2+3159 11 



20 



15 



10 




4500 



5000 5500 6000 

Rest Wavelength (A) 



6500 



Figure 7. Subtraction of Fe II template from the observed spectrum of RX J1334. 2+3759. The vertical scale is the relative flux. The 
top spectrum is the redshifted observed spectrum. The bottom one is the Fc II model for RX J1334. 2+3759. The middle spectrum is the 
Fe II subtracted spectrum of RX J1334. 2+3759. Note that the Fe II corrected spectrum has been shifted downward for visibility. 



in the density. The inner regions of the BLR in the NLS1 
galaxies can be expected to be fully ionized in hydrogen due 
to higher ionizing power of the steeper ionizing continua, 
and the temperature there can be higher than in the inner 
BLR of the Seyfert 1 galaxies. These kind of physical con- 
ditions in the inner BLR of the NLSls will not, however, 
produce an appreciable number of H/3 photons, while the 
outer regions of the BLR will have favorable physical con- 
ditions to generate a copious number of H/3 photons. UV 
emission lines, which probe the inner BLR, are expected to 
be broader than the H/3 line. Thus, an enhancement of the 
density due to outflows can explain (a) the lower width of 
H/3 and Ha in the NLSls compared to Seyfert 1 galaxies, 
and (b) comparable width of the permitted lines in the UV 
spectra of NLS1 galaxies compared to Seyfert 1 galaxies. It 
should be noted that a relatively high density (Jj 10 11 cm ) 
BLR has been inferred from the UV spectrum of I Zw 1 
- a prototype NLS1 object (Laor et al. 1997b). There is 
also observational evidence for outflows in the NLS1 galax- 
ies e.g., a weak UV absorption system with a line of sight 
outflow velocity of ~ 1870 km s _1 has been detected from 
I Zw 1 (Laor et al. 1997b ). Outflows in three other NLS1 
galaxies, viz., Akn 564, WPVS007 and RX J0134.2-4258, 
have also been observed by the presence of UV absorption 
lines in their spectra (Goodrich 2000). In addition, Leighly 
ct al. (1997 ) have found evidence for relativistic outflows in 



the ASCA spectra of 3 NLS1 galaxies. Higher accretion rates 
that would likely produce stronger outflows thus pushing the 
BLR further radially outward thereby resulting in narrower 
Balmer lines in the NLS1 objects, also result in steeper soft 
X-ray continua. Thus, a variation in the higher accretion 
rate can result in the observed anti-correlation between the 
slope of the soft X-ray continua and the width of the H/3 
line in the NLS1 objects. 

Another advantage of the above picture is that the 
strong outflows can result in shocks which may be respon- 
sible for the Fe II emission observed in the NLS1 galaxies. 
Collin & Joly (2000) have shown that the standard photoion- 
ization models with any set of parameters cannot produce 
the observed strength of the Fe II emission in the NLS1 
galaxies, and have suggested that non-radiative heating, for 
instance, shock heating for the production of Fe II emission. 



5 CONCLUSIONS 

(i) RX J1334. 2+3759, a newly discovered narrow-line 
quasar, is highly luminous in soft X-rays (Lx ~ 2.2 x 
10 44 erg s _1 in the energy band of 0.1-2.0 keV). Soft X-ray 
emission from RX J1334. 2+3759 is very steep (Fx ~ 3.8) 
and is highly variable. The most extreme variable event has 

erg s~ 2 



^£ = (1.95 + 1.02) x 10 42 
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Figure 8. Decomposition of observed H/3 profile into two Gaussian components. Vertical scale represents relative flux. 
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Figure 9. Decomposition of Ho profile into two Gaussian profiles. Vertical scale represents relative flux. 
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(ii) The optical spectrum of RX J1334. 2+3759 is typical of 
NLS1 galaxies. Decomposition of Balmer H/3 and Ha lines 
has revealed the presence of narrow (FWHM~880 kms" 1 ) 
and broad components (FWHM ~ 2850 km s _1 ) in them. 

(iii) The ratio of [O IIIJA5007 line flux and the flux in the 
narrow component of H/3 line is ~ 0.30 which is very dif- 
ferent from the value of ~ 10 found in Seyfert 2 galaxies 
indicating that the NLR in RX J1334.2+3759 is different 
from the NLR in the normal Seyfert galaxies. 

(iv) A possible explanation for the observed properties of 
the NLR and BLR of RX J1334.2+3759 and other NLS1 
galaxies, as well as for the well known anti-correlation be- 
tween soft X-ray slope and H/3 width found in them, has 
been suggested in terms of density enhancements resulting 
from outflows due to super Eddington accretion rates. 
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